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Tools  for  Time-Frequency  Analysis 


Frames  of  windowed  Fourier  transforms  and  wavelets;  work  in  collaboration  with  Radu  Balan,  a 
graduate  student  supported  by  this  grant. 

Frames  provide  redundant  descriptions  of  vectors  in  a  Hilbert  space;  that  is,  the  sequences  (</,  (pk>)keK. 

where  the  (tpiJksK  constitute  a  frame  in  ^  and  /  roams  over  all  of  do  not  fill  all  of  ^^(K)  because  of 
linear  dependencies  between  the  (pk  reflected  by  the  inner  products  </,  (pk>-  The  range  of  the  frame,  i.e. 

the  subspace  of  ^^(K)  given  by  {(</,  (Pk>)keK;  characterizes  the  frame  up  to  an  isomorphism:  if 

two  frames  (9k)keK  and  (\|/k)keK  have  the  same  range,  then  there  exists  a  banded  operator  A  with  banded 
inverse  such  that  \|/k  =  Atpk-  This  and  more  detailed  observations  led  to  [1],  a  paper  studying  the  geometric 
structure  of  frames. 

On  the  other  hand,  if  the  ranges  of  two  frames  are  orthogonal  subspaces  of  ^^(K),  then  we  call  the  frames 
themselves  orthogonal.  In  this  case  it  is  possible  to  use  them  for  multiplying.  This  means  that  if  /,  g  are 

two  (unrelated)  functions  in  we  can  recover  both  f  and  g  from  the  single  sequence  (ak)keK.  with  ttk  = 
</.  9k>  +  <8,  Vk>-  This  construction  can  be  extended  to  the  case  where  the  ranges  of  the  two  frames 
make  a  non-zero  angle  [2]. 

In  a  dual  approach,  we  are  also  investigating  how  frames  can  be  used  for  multiple  descriptions.  In  this 
case,  a  redundant  description  of  a  signal  gets  split  into  different  channels,  each  of  which  carries  only 
partial  information.  The  idea  is  then  to  design  the  scheme  optionally  so  that  even  if  some  of  the  channels 
are  interrupted,  the  reconstruction  from  the  remarking  partial  information  is  best  possible.  This  is  work  in 
progress,  with  Vinay  Vaishampayan. 

In  addition,  Radu  Balan  also  proved  several  stability  results.  For  instance,  he  proved  that  for  any  wavelet 
Riesz  basis,  one  can  stretch  or  shrink  the  translation  parameter  slightly  and  still  retain  a  Riesz  basis  [3]. 

Time-frequency  localization: 

The  squeezed  wavelet  transform  used  in  the  joint  work  with  S.  Maes  [4]  is  a  special  case  of  a 
reassignment  method.  With  Patrick  Flandrin  and  Eric  Chassande-Mottin,  the  PI  has  studied  other 
reallocation  schemes.  Here  one  starts  from  a  time-frequency  decomposition,  in  which  one  tries  to  “undo” 
the  blurring  effect  of  the  time-frequency  tool  so  as  to  obtain  more  sharply  focused  time-frequency 
representations.  We  have  succeeded  in  finding  some  forms  that  reallocate  in  both  the  time  and  frequency 
variables,  while  still  giving  exact  reconstruction. 

Moreover,  we  have  studied  a  “differential  reassignment  scheme”  that  lets  a  signal-dependent  vector  flow 
act  on  the  Wigner  transform;  under  this  flow,  the  transform  concentrates  into  attracting  basins,  which 
makes  it  possible  to  use  this  as  a  tool  to  decompose  signals  into  different  components  [5].  In  a  different 
approach,  we  are  exploring,  with  Fabrice  Planchon,  the  mathematical  aspects  of  smoothing  the  Wigner- 
Ville  distribution  with  Gaussians  that  are  locally  oriented  and  stretched  in  a  signal-adaptive  way,  given  by 
geometric  properties  of  the  Wigner  distribution  itself.  This  results  in  a  positive  time-frequency 
representation  without  the  interference  terms  that  characterize  Wigner  distributions,  but  that  is 
nevertheless  sharper  than  most  smoothed  Wigner-Ville  distributions. 

Lossless  compression: 

In  some  applications,  the  data  that  one  wants  to  compress  are  already  in  the  form  of  integers  (digitized 
grey  scale  images,  for  instance).  With  A.R.  Calderbank,  W.  Sweldens,  and  B.L.  Yeo,  the  PI  developed 
schemes  that  adapt  wavelet  transforms  into  schemes  that  map  integers  to  integers.  These  turn  out  to  be 
generalizations  of  schemes  that  exist  in  the  literature.  We  can  use  our  method  starting  from  any  wavelet 


filter  bank,  and  for  some  of  these  filter  banks  we  obtain  better  results  for  subband-based  lossless  image 
compression  than  previously  known  [6]. 

Recently,  B.L.  Yeo  managed  to  push  these  results  even  further,  and  our  approach  now  leads  to  better 
results  than  any  other  lossless  scheme,  whether  or  not  it  uses  subband  filtering  [7]. 

Smoothness  for  unequally  spaced  subdivision: 

With  Igor  Guskov  (Ph.D.  student  partially  supported  by  this  grant)  and  Wim  Sweldens  (Bell 
Laboratories)  we  studied  the  smoothness  for  subdivision  schemes  with  non-equally  spaced  points. 
Subdivision  schemes  underlie  every  construction  of  wavelet  bases  linked  with  fast  algorithms.  The 
smoothness  generated  by  the  subdivision  scheme  is  the  key  to  the  smoothness  of  the  wavelets  as  well  as 
to  their  potential  for  data  compression,  or  for  the  sparse  yet  accurate  representation  of  multiscale 
phenomena.  Although  we  have  a  very  thorough  and  detailed  understanding  of  these  smoothness  issues  in 
the  case  where  the  wavelet  centers  are  placed  on  regularly  spaced  lattices,  very  little  is  know  for  the  case 
where  the  wavelet  centers  are  irregularly  spaced.  A  better  understanding  of  these  issues  is  crucial  for 
many  applications,  including  compression,  visualization  or  editing  of  signals  or  surfaces  given  by 
irregularly  sampled  data.  We  have  shown  that  although  neither  Fourier  nor  spectral  methods  apply  in  the 
irregularly  spaced  case,  we  can  nevertheless  still  use  “commutation  formulas,”  which  are  constraints  on 
the  subdivision  filters,  to  obtain  smooth  functions  [8]. 
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Interactions: 

Together  with  Stuart  Schwartz  (Electrical  Engineering)  and  Rene  Carmona  (Civil  Engineering),  and  now 
also  joined  by  Michael  Orchard  (Electrical  Engineering),  the  PI  organized  a  weekly  brown  bag  seminar  on 
time-frequency  and  time-scale  methods,  featuring  inside  as  well  as  outside  speakers.  Participants  are  faculty, 
postdocs,  and  graduate  students  from  Electrical  Engineering,  Civil  Engineering,  Mathematics,  Applied  and 
Computational  Mathematics,  and  Computer  Science.  These  meetings  were  very  lively,  and  usually  extended 
well  into  the  afternoon  with  interesting  discussions.  These  often  led  to  the  formulation  of  new  problems  to 
look  at;  they  also  gave  rise  to  new  collaborations.  We  are  also  considering  organizing  one  or  two  one-day 
meetings  of  the  same  nature  with  a  wider  circle  of  time  frequency  researchers  in  the  Northeast. 
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